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Abstract 

We consider critical curves — conformally invariant curves that 
appear at critical points of two-dimensional statistical mechanical sys- 
tems. We show how to describe these curves in terms of the Coulomb 
gas formalism of conformal field theory (CFT). We also provide links 
between this description and the stochastic (Schramm-) Loewner evo- 
lution (SLE). The connection appears in the long-time limit of stochas- 
tic evolution of various SLE observables related to CFT primary fields. 
We show how the multifractal spectrum of harmonic measure and other 
fractal characteristics of critical curves can be obtained. 
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Figure 1: A critical Ising model. Black and white represent ±1 
spins. The cluster boundaries furnish an example of critical curves (after 
http: / / www.ibiblio.org/ e-notes /Perc /contents. htm) . 



1 Introduction 

The continuous limit of two-dimensional (2D) critical statistical systems 
exhibits conformal invariance which proved to be a very useful and compu- 
tationally powerful concept [H [2j [3] . Because of the conformal invariance a 
system with a boundary can be studied in any of the topologically equivalent 
domains. We will consider critical systems in domains with the topology of 
an annulus which for the ease of computations can be mapped by a confor- 
mal transformation onto the upper half plane with a puncture. 

Critical 2D systems can be described in terms of fluctuating curves 
(loops, if they are closed) which can be understood as domain walls or bound- 
aries of clusters (e.g. Fig. [I]). The exact definition depends on the system 
in question and typical examples include the boundaries of the Fortuin- 
Kasteleyn clusters in the g-state Potts model, the loops in the 0(n) model, 
and the cluster boundaries in a critical percolation system. We use the term 
critical curves to describe any of these objects. 

Critical curves are drawn from a statistical thermal ensemble and because 
of the scale invariance at criticality they are self-similar random fractals. 
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Figure 2: A critical curve is conditioned to start from 0, pass through a 
puncture and end at L on the boundary. Another critical curve is condi- 
tioned to end at a puncture in the bulk. 

A large scope of their geometrical properties can be investigated. As an 
example, one can study the scaling of the electric field near a critical curve 
which is assumed to be the boundary of a conducting and charged cluster. 
Averaged over the fluctuations of the curve (we denote such averages by 
-<...>- ) , the moments of the electric field E scale with the distance r to 
the curve as -<E h (r)>- ~ C(h)r A ^ with a non-trivial multifractal exponent 
A(h) first computed in [H[5]. On the one hand, this exponent is related to 
the multifractal spectrum of the harmonic measure of the critical curve (and, 
in particular, determines its fractal dimension). On the other hand, A(/i) 
happens to be the gravitationally dressed dimension h introduced in [7], thus 
forming an interesting link to quantum gravity. The determination of the h- 
dependent prefactor C{h) is also interesting and, to our knowledge, has not 
been done yet. One can also study more complicated correlators of moments 
of the electric field measured at different points: -<E hl (zi)E h2 (22) • ■ •>- (we 
use complex coordinates z to denote points in the plane). Such objects probe 
local properties of critical curves. 

In this paper, which extends our results published in Ref. [6], we consider 
both the CFT and SLE approaches to stochastic geometry of critical curves 
and elaborate on the relation between the two. We also obtain the results 
of [H [5] for multifractal exponents using CFT without reference to quantum 
gravity. 
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2 Creation of critical curves and conformal invari- 
ance 

As the statistical system fluctuates critical curves come and go in a random 
fashion and this makes them difficult to study. We can, however, choose a 
point in the bulk or on the boundary and remove from the partition function 
all realizations except for those in which there is a curve passing through this 
point. The physical meaning of such conditioning is transparent. Consider 
a critical Ising model in the upper half plane and define critical curves to 
be the boundaries of clusters of spins of the same sign. As a boundary 
condition, we demand that all spins on the negative real axis are +1, and 
all those on the positive real axis are —1. This condition guarantees that in 
all realizations there will be a critical curve growing from the origin. 

This is just one of many imaginable ways of conditioning the system but 
it illustrates the point: the existence of a curve growing from a point on the 
boundary is ensured by a change of the boundary condition at this point. 
To condition the curve to pass through a bulk point one can insert there a 
puncture, thus effectively turning it into a boundary point (Fig. [2]). 

Critical curves are conformally invariant [H [2] (and are sometimes also 
called conformally invariant curves) in the sense depicted in Fig. [3j Define 
a system in two different but topologically equivalent domains A and B 
and introduce a conformal transformation f(z) which maps A onto B. We 
consider the ensemble of loops in both domains and also the ensemble C 
of images of loops produced by the conformal map f(z). The lattice in C 
is warped by the conformal transformation which makes it different from 
the system in B. Nevertheless, it may happen that the loops in B and C 
are statistically identical to each other, that is, they are different random 
realizations drawn from the same statistical ensemble. In this case we say 
that the ensemble of loops possesses conformal invariance. 

The continuous limit of critical systems with conformal invariance is de- 
scribed by conformal field theory (CFT) which is characterized by a central 
charge c. In this paper we consider only theories with c ^ 1. For — 2 ^ c ^ 1 
they can be obtained as a continuous limit of the lattice 0(n) model. 

The central charge itself does not completely determine the critical curves: 
their behavior crucially depends on boundary conditions. A CFT with c ^ 1 
can be represented as a theory of a Bose field with the Neumann or the 
Dirichlet boundary conditions. From the point of view of the 0{n) model, 
different boundary conditions correspond to dense and dilute phases in which 
the behavior of critical curves is different. An impressive progress in under- 
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z -> f(z) 



Figure 3: A and B are critical systems defined in different domains. The 
system C is the image of A under the conformal transformation f(z) which 
maps domain A onto domain B. Conformal invariance is the statistical 
identity of B and C. 




Figure 4: Schematic behavior of dilute (left) and dense (right) loops. The 
external perimeter of the dense curve is a dilute curve. 

standing stochastic geometry of critical curves was achieved in the stochastic 
Loewner evolution (also known as Schramm-Loewner evolution, or SLE) [8]- 
|19j . It can be deduced that curves in the theory with the Dirichlet boundary 
condition (dilute phase) are almost surely simple. This means that the prob- 
ability for them to touch themselves or each other vanishes. Those with the 
Neumann boundary condition for c —2 (dense phase) almost surely touch 
themselves and each other (Fig. For c < —2 the curves are plane-filling. 

There exists a duality relation between dense and dilute phases |20j . 
In particular, the external perimeter of a dense curve is a dilute curve [5]. 
Thus, studying the dilute phase also answers some questions in the dense 
one. In CFT this duality is expressed by the fact that both phases can have 
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the same central charge and the difference lies in the boundary conditions. 



3 Boundary operators and stochastic geometry 

Consider a system in the dilute phase occupying the upper half plane and 
conditioned to contain a critical curve 7 connecting the points and L on 
the real axis (see Fig. [2]). As was mentioned above, this will be the case if the 
boundary condition changes at these two points. Such a change is produced 
by the insertion of a special boundary condition changing operator which we 
can call a curve-creating operator. In statistical models with c ^ 1 such as 
the 0{n) model it is known [21] and we will discuss it later. The correlator 
of two such operators 

(B(0)B(L)) m = £M 

is determined by the partition function Z(0, L) with boundary conditions 
changing at and L. Here Z is the partition function of the system with 
uniform boundary conditions. The subscript of the correlator always refers 
to the domain of definition. We note in passing that by fusing together sev- 
eral such ,6's we obtain other boundary condition changing operators whose 
insertion produces multiple curves growing from a point on the boundary. 

This correlation function can be interpreted in terms of stochastic geom- 
etry of the curve. The interpretation involves a two-step averaging (similar 
arguments can be identified in [lOl [22] )• I* 1 the first step we pick a par- 
ticular realization of the curve 7. Then 7 is the boundary separating two 
independent systems — the interior and the exterior of 7. In both these 
systems we can sum over the microscopic degrees of freedom to obtain the 
partition functions Zi nt and Z^ xt , respectively. These are stochastic objects 
that depend on the fluctuating geometry of 7. 

In the second step we average over the ensemble of curves of 7. We thus 
obtain 

Z(0,L) = ^Z^Z^y , 

where -<. . .>- stands for averaging over the shape of 7. 

We further consider the insertion of an additional field O(z) on the 
boundary and a distant field \P(oo) such as to make non-zero the correlation 
function 

(B(0)O(z)B(L)t>(oo)) m (1) 
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The exact form of \l/(oo) is not important, for example we can choose it 
as 0(oo). If we are only interested in the z-dependence of this correlation 
function in the limit \z\ <C \L\, we can fuse together the distant fields: 
jB(L)^(oo) :— > ^(oo) and therefore consider a three-point function 



In the two-step averaging procedure we rewrite dU) as an average over the 
fluctuating geometry of 7: ^(C(z)*(oo)) 7 Zj ) nt Z^ x V , where (C(z)*(oo)) 7 
stands for the correlation function in the exterior of 7. In the limit \z\ <C 
\L\ this correlation function is statistically independent from the other two 
factors, and we are left with the scaling relation: 



The two-step averaging can be extended to multiple operators. If we 
choose O(z) a primary field, this relation yields the statistics of harmonic 
measure of critical curves (see Sec. [6]). 

4 Loop models and the Bose field 
4.1 Dense phase 

The relation between critical curves and operators of a CFT with a boundary 
is most transparent in the representation by a Bose field tp(z, z) [3ll20 l l23 l l2l]. 
This representation is commonly known as the Coulomb gas method. As 
in the case of any field theory, its derivation from the lattice model is not 
rigorous but is a posteriori justified by the results. 

A large variety of 2D lattice statistical models (percolation, g-state Potts, 
0(n), XY, the solid-on-solid and other similar models) can be defined in 
terms of a system of random curves, usually closed loops unless special 
boundary conditions are imposed. Each loop has a fixed statistical weight 
in the ensemble. Mappings between specific lattice models and loop models 
are described in detail in many reviews |20[ [25l l26| [27] . Here we consider the 
0(n) model as a representative example. It covers a large range of central 
charges (— 2 ^ c ^ 1) and its dilute and dense phases correspond to the 
Dirichlet and Neumann boundary conditions of CFT. The model can be 
defined on the honeycomb lattice on a cylinder (such as a carbon nanotube) 
in terms of closed loops (Fig. [5]): 



(fl(O)0(z)tf(oo)>H 



(jB(0)O(z)jB(L)*(oo)) h ~ ^(0(z)*(oo)> 7 ^ 




(3) 



loops 
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Here a; is a parameter related to the temperature, L is the total length of 
all loops and N is the number of loops. The parameter n is not necessarily 
an integer. 

The 0(n) model has a critical point at 

x = x c (n) = — 

for all n in the range — 2 ^ n ^ 2. At the critical point the typical length 
of a loop diverges but the loops are dilute in the sense that the fraction of 
the vertices visited by them is zero. For x > x c (n) the typical length of the 
loops still diverges so that they are still critical but visit a finite fraction 
of the sites. This is the dense phase of the model. In the continuous limit 
dilute loops are simple and dense ones can touch themselves and each other. 
Finally, at zero temperature (x = oo) the loops go through every point on 
the lattice. This is the fully packed phase. In the continuous limit the loops 
become plane filling. 

The continuous description of the critical behavior of the 0(n) model 
by a local field theory is best understood in the dense phase. We need to 
have a description of the system in terms of local weights on the lattice and 
this is done as follows. One randomly orients the loops and assigns a weight 
to each orientation. In order to reproduce the partition function ([3]) the 
sum of the weights for two orientations of every loop should be n. This is 
achieved by giving a local complex weight e ±ieo7r / 6 to each lattice site where 
an oriented loop makes a right (left) turn. The weight of an oriented closed 
loop is the product of weights on all sites that it visits. For a closed loop 
which does not wrap around the cylinder this product equals e ±teo7T since 
the difference between the numbers of right and left turns is ±6. The sum 
over the two orientations gives the real weight n for an un-oriented loop if 
we choose 

n = 2 cos 7reo- 

The range of —2 ^ n ^ 2 is covered once by ^ eo ^ 1. To describe both 
the dilute and the dense phases, however, we will need to allow for a wider 
range —1 ^ eo ^ 1, with positive eo for the dense phase and negative eo for 
the dilute phase. 

The situation is different for the loops which wrap around the cylinder. 
For them the difference between the numbers of right and left turns is 
and they are wrongly counted with weight 2 instead of n. To correct this, 
one defines for each configuration of oriented loops a real height variable H 
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Figure 5: If a loop wraps around the cylinder, the difference between the 
numbers of left and right turns is 0, otherwise it is ±6. On a flat domain of 
the same topology — annulus, — this difference is ±6 for all loops. 



which resides on the dual lattice and takes discrete values conventionally 
chosen to be multiples of ir. We start at some reference point where we 
set H = 0, and then every time we cross an oriented loop, we change H 
by ±7r depending on whether we cross the loop from its left to its right or 
vice versa. Then one introduces an additional factor e ie o(H{+oo)~H(~oo)) - m ^ 
each term of the partition function, where ±oo stand for two chosen points 
on the two boundaries of the system (the bases of the cylinder). This factor 
is 1 if there are no wrapping loops, so it does not affect the weights of non- 
wrapping loops. But since each wrapping loop changes H(+oo) — H(—oo) 
by ±7r, this factor turns into the proper weight n for each wrapping loop. 
Thus, we represent the partition function ([3]) in the form 



We note that the insertion of the fields at the boundaries is brought about 
not only by the topology of the domain: in a flat annulus, although its 
topology is the same as cylinder's, all loops are counted with the weight n 
as it is, so that no extra factors are needed (Fig. [5]). 

When we pass to the continuous limit, the lattice height function H is 



.cyl 
O(n) 





(4) 



oriented loops 



10 



coarse-grained and becomes a fluctuating compactified scalar field h(z,z): 

h~h + 2n. 

Upon coarse-graining the 0(n) loops on a cylinder become level lines of the 
field h. In the dense phase this field is believed to be Gaussian, that is, its 
action contains the term (g/Air) J d 2 x (V/i) 2 , where the fluctuation strength 
parameter g is yet to be determined. The boundary terms of have to 
be added in any geometry in the form of an oriented boundary integral 
(ieo/27r) Jg D dl Kh where K is the geodesic curvature of the boundary. This 
is consistent with (HJ). 

A similar situation occurs if the critical system lies on a surface with 
curvature which microscopically can be viewed as existence of defects on 
the honeycomb lattice (pentagons and heptagons correspond to positive and 
negative curvature, correspondingly). The correct weight for a loop that 
surrounds a region of non-zero Gaussian curvature R is obtained if we include 
in the action the so-called background charge term (ieo/87r) J D d 2 xRh. 

The determination of the coupling constant g can be done either by 
comparison with the exact solutions of a six-vertex model related to the 
0(n) model or by an elegant argument of |28^ I29j. Another necessary term 
in the action is the locking potential of the form A j d 2 x V(h) which favors 
the discrete values h = 0, 7r. It must be, therefore, a 7r-periodic function 
of h, the most general form of it being V = X^ez i^o vie 2llh . We treat the 
locking potential perturbatively. In the zeroth order (A = 0), each term of 
V is a vertex operator of dimension j20| xi = (2/g)l(l — eo). An operator 
with x\ < 2 is relevant. For the perturbative approach to be applicable, the 
locking potential should contain no relevant terms. Furthermore, the most 
relevant term should be marginal. This fixes g. For < eo < 1 the most 
relevant term is I = 1 and this gives the relation 

eo = l — g, n = —2cosirg. (5) 

Note that < g < 1, as it is expected in the dense phase. 

This description is not valid in the dilute phase. There H does not renor- 
malize into a Gaussian field. A manifestation of this is that it is impossible 
to make the locking potential marginal for 1 ^ g ^ 2. Naively, this requires 
taking — 1 < eo < 0, but then one needs to pick I = — 1 term as the most 
relevant and it still gives us g = 1 + eo = 1 — ] eo ] < 1. 

The point g = 1 separates the phases and is somewhat special: there the 
/ = 1 and I = — 1 terms in the locking potential have the same dimension 
and we have to keep them both. The relation between Gaussian fields and 
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critical curves has been made rigorous in |30j which in our language describes 
the case n = 2, g = 1. 

We now introduce the parametrization 

4 

9 = -■ 

K 

The range < g ^ 2 where the 0(n) model exhibits critical behavior corre- 
sponds to 2 ^ k < oo. Smaller values < k < 2 correspond to multi-critical 
points in the 0(n) model but also make sense for other critical systems. 
Notice that k < 4 and k > 4 describe the dilute and the dense phases cor- 
respondingly, while k = 4 gives the point g = 1 separating the two phases. 
Readers familiar with SLE will recognize that this is similar to the SLE 
phases [31]. In fact, k can be identified with the SLE parameter. For k ^ 4 
SLE curves are simple and for 4 < k < 8 they have double points. 

In the literature it is customary to rescale the field <p = ^/2gh and fix the 
coupling constant g to be 1/2 at the expense of varying the compactification 
radius of ip: 



R =Vk (6) 

This is the normalization that we adopt from now on. Up to the locking 
potential one then writes the partition function as 



" v / 2 r "" 

27T J dD 



Z oin) ~ / IV cos $ dlKcp e -/(^) 2 fr, ( 7 ) 



where the relation ([5]) reada3 



Jk 2 
2a = \ > 0. 

2 \ K 



The marginal locking potential now reads A J d 2 xe % ^ a+Lp , where we used 
the notation 



Jk 2 

q+ = — , a_ = -=. (8) 

2 V K 

The partition function written in this way remains valid for any flat geometry 
with any number of boundaries. 



1 The notation differs from [6], where ao < in the dense phase. As a consequence, in 
the Kac notation ^v jS of curve-creating operators r <-> s. 
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Note that a± satisfy simple relations a± = a® ± \J a§ + 1. The locking 
potential is treated as a perturbation: in each order of the perturbation 
theory a number of screening charges j d inserted into all cor- 

relation functions so as to satisfy the neutrality condition [32]. Note that 
the second screening charge f d 2 x e l ^ a ~ v , although also of zero dimension, 
is not found by the argument on the lattice — it is incompatible with the 
compactification radius. 

If the domain's Gaussian curvature R is not zero (as it happens when 
heptagons or pentagons are inserted as defects of the lattice) an additional 
term i v ^"° J d 2 xRtp is needed. 

Due to the symmetry (p — > —99 we can rewrite the partition function 
as a local field theory Zo( n ) ~ f T)ipe~ s , but with a complex action 5 = 
A + A f D d 2 x e iv/ 2 Q +^ where 

1 f \J~2oi. f 
A = — / d 2 x (V(p) 2 + i / dl K<p, dense phase, «o > 0. (9) 

Jd 2n JdD 

Since this action is local we can use the usual methods of conformal field 
theory. Such an action is known to describe a field theory with central 
charge 

2 (« - 4) 2 
c = 1 - 24ag = 1 - 3^ '—. 

However, the saddle point of the action does not lie in the real fields and 
therefore the Bose field in the OPEs is not real so that it is no longer the 
coarse-grained height function. 

4.1.1 Conformal transformations and boundary conditions 

As an illustration, we consider the upper half plane (with punctures). In 
this case the curvature of the boundary is zero except at the point at infinity. 
By conformal invariance, the theory can be mapped onto any domain with 
a curved smooth boundary. 

Since the field is Gaussian it consists of the holomorphic and the anti- 
holomorphic parts: 

ip(z,z) = 4>(z) + 4>(z), 
whose correlation functions are 
(0(z)^ / )) = -log(z-z / ), $(zW)) = -l0g(z-it), (0(s)0(z)> = O. (10) 
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We will also use the dual boson field <p related to ip through Cauchy-Riemann 
conditions: 

(f(z,z) = — i4>(z) + i(f)(z). (11) 

Varying the action A with respect to the metric we find the components of 
the stress-energy tensor: 

T(z) = ~(d<P) 2 + iV2a o d 2 0, f(z) = ~(34>) 2 + iV2a o d 2 0, (12) 

where d = d/dz,d = d/dz, and the products are normal ordered. 

When the critical system is defined in the upper half plane the con- 
formally invariant boundary condition is T{x) = T(x) for all s £ 1 |21] . 
For the stress-energy tensor (|12|) this condition readily yields 4>(x) = (fi(x), 
therefore (p = 4>{z) + 4>(z) and (p = —i(fi(z) + i(fr(z) which is equivalent to the 
Neumann boundary condition: 

^1^=0 = °- 

The Neumann condition implies that (p is a scalar which agrees with the 
meaning of the Bose field as height function. 

The components (|12|) determine the transformation properties of the 
fields under a conformal transformation z — > f(z): 

^)^ ( />(/(z)) + iv / 2a log/ / (z), (13) 

which yields 

<p(z, z) <p(f(z),W)) + *2v^a log \f'(z)\. 

This transformation is merely the transformation of the classical configura- 
tion of ip (the saddle point of (0). Notice that although <p is not real, its 
real part does not transform so that the level lines 

Reip(z, z) = const. 

are conformally invariant fluctuating curves. This suggests to identify them 
with the continuous limit of critical curves i.e. the level lines of the coarse- 
grained height function. 
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4.1.2 Electric and magnetic operators 

In the following we will see that critical curves in both dense and dilute 
phase are created by spinless vertex operators: 

0^ e ' m \z,z) = ^^ e H>{z,z) e -s/2nup{z,z) ^ Q4) 

where we introduced "electric" and "magnetic" charges e and m. The 
boundary term in the action is the insertion of an electric background 
charge e = 2ao at a puncture in the upper half plane. 

A vertex operator is a primary field and can be written as a product 
of the holomorphic and the antiholomorphic components 0^ e ' m \z, z) = 
V a (z)V a (z), which are assumed to be normal ordered: 

V a (z) = e i ^ a< l > ( z ) i V & (z) = e * v/ 2 |5 ^( 5 ). 

Here the holomorphic and antiholomorphic charges are 

a = e + m, a = e — m. (15) 

In the dense phase the Neumann boundary condition allows to write 

0^ e,m \z,z) = e i v / 2(e+m)</)(z) e iv / 2(e-m) 9 i(5)_ Qg\ 

This can be interpreted as gluing together the holomorphic and the an- 
tiholomorphic sectors which is best understood in terms of image charges 
[2TJ [231 121] • The operator (fT6|) is regarded as a product of two holomorphic 
operators: one in the upper half plane and the other at the image point 
in the lower half plane. In the image the electric charge remains the same 
while the magnetic charge changes sign. 

The holomorphic and the antiholomorphic weights of the vertex opera- 
tors are found by applying the stress-energy tensor (fT2|) to them: 

h a = a (a — 2ao), ha = a(a — 2chq). 

A given weight h corresponds to two charges: ao ± \J (Xq + h. 
In terms of the electric and magnetic charges the weights read 

h(e, m) = (e + m)(e + to — 2ao), h(e, m) = (e — m)(e — to — 2olq). (17) 

A vertex operator is spinless (meaning that h = h) if either a = a or 
a = 2ao — a. In the first case the operator is purely electric (to = 0) and in 
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the second case it can have an arbitrary magnetic charge m but the electric 
charge should be e = ctQ. 

We will also find it convenient to denote the vertex operators by their 
weights. Thus Oh(z, z) is a spinless vertex operator of weight h in the bulk 
and Oh(x) is a boundary vertex operator of weight h. 

The propagators (fT0|) lead to fusion by the addition of charges: 

lim V ai ( Zl )V a2 (z 2 ) = Oi - z 2 ) 2aia2 V ai+a2 (z 1 ) + .... (18) 

Z\ — >Z2 

In the literature it is customary to label holomorphic charges and weights 
by two numbers r, s according to 

a r ,s = a - ^(™+ + sa_) = -(1 - r)a + + -(1 - s)a_, 

. . (rn - 4s) 2 - (k - 4) 2 
"r, s = a r , s a r , s - 2a = 77 • (19) 

lOK 

The primary field of weight /i r ,s is denoted tpr,s- We will use it as a shorthand 
notation without imposing any restrictions on r or s. 



4.1.3 Curve-creating operators in the dense phase 

In this section we determine the vertex operators which create critical curves. 

On the lattice, the appearance of n critical curves emanating from a 
small region in the bulk results in the change in the height function H by 
±7rra as one makes a full circle around this region. This assumes that all 
of the curves are oriented the same way — inwardly or outwardly, — for 
otherwise they would be able to reconnect with each other. In the continuous 
limit H — ► Re ip this represents a vortex configuration of the field (p. As the 
propagators (|10p show, a vortex is produced by the insertion of a magnetic 
charge: 

0( e -°) (z, z)Q^ m ) ( z ', z') = ( ^—4^ 2e "V e + m ) (z')0 (e " m) (z') + ... 

\Z — Z ' 

— e 4iem axg(z- z') Q(e+m) ^ z '^Q( e - m ) + . . . 

This means that for z — > z' the value of the field changes by ±4\/2vrm as 
one goes around z' in full circle. Hence, a discontinuity line arises (Fig. [6]). 
If this vortex corresponds to a star of n critical curves joining at the point 
z' , the change in ip should be equal to ±mrlZ. The compactification of ip 
makes the discontinuity line invisible for an even number of curves. For an 
odd number of curves, a discontinuity line with the jump ttTZ remains. One 
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can view it as a critical curve whose form is conditioned to be fixed so that 
it is removed from the ensemble of random curves. 

We then find the magnetic charge of the bulk operator which creates n 
curves as 

y/2 _ n n 

where we used the value (jSJ) of the compactification radius and the definition 
([5]) of a_. In order to be spinless (otherwise the operator would transform 
under rotations, giving non-trivial dependence on the winding number of 
curves) , the bulk curve-creating operator should also have an electric charge 
ceo thus acquiring the form 

V ao+m (z)V a °- m (z). 

We have therefore found the possible holomorphic charges of the bulk curve- 
creating operator to be 

n 

" = "o =F -7OL- = a 0i±n / 2 

in the notation Q19|). These charges correspond to the same weight 

4n 2 - (« - 4) 2 



h 



0,71/2 



16k 



The curve-creating operator itself is then denoted ^o,n/2- I n particular, a 
single critical curve going through a point z (which is the same as n = 2 
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critical curves meeting at the point z) is created by the operator -00, 1 with 
the holomorphic weight /io,i = (8 — re)/16. Note that this weight determines 
the fractal dimension of a critical curve as 2 — 2/&o,i [33J. 

The creation of curves on the boundary is considered in the same way. If 
n curves exit a boundary point the value of the height function changes 
by ±nirlZ as one goes around this point in a semi-circle. On the real 
axis [e ±iV ^ a ^ x \ ^)(y)} = ±^ita9{x - y ) e ±^<t>{*) where 9(x) is the step- 
function. This shows that the field value (p on the boundary changes by 
in crossing the inserted vertex operator. Moreover, in presence of 
the Neumann boundary condition p(x) = 2<p(x) on the real axis. The n 
curves are created, therefore, by the insertion of a boundary operator of the 
form e ± i v / 2"i,n+i<)! , _ The weight depends on the sign in the exponent. The 
lower weight corresponding to plus sign so that this operator is more rel- 
evant. This is then the curve-creating operator on the boundary denoted 
ipi n+i- There also exist operators which condition a curve not to touch a 
certain part of the boundary [331 [Ml l35| 155]. 

4.2 Dilute phase a < 

Dilute critical curves occur, for example, in the Oin) model at the critical 
temperature. They are known to be described by g > 1 (k < 4). The 
central charge and the conformal weights of all curve-creating operators 
are obtained by continuation of the corresponding formulas from the dense 
phase (« > 4). The arguments used to derive the action of field theory from 
the lattice in the dense phase cannot be straightforwardly extended to the 
dilute phase. Nevertheless, the description of the dilute phase is available 
via the electric-magnetic duality transformation e <-> m also called T-duality 
in the literature (see [5J [20l [371 [38] ) ■ Every correlation function of primary 
fields — and in particular curve-creating operators, — can be represented 
as the insertion of electric and magnetic charges into the partition function: 
Z[g,e,m]. The electric-magnetic duality means Z[g,e,m] = Z[l/g,m,e]. 

The interchange between the electric and magnetic background charge 
means the following change of the boundary term in the action ([9]): 

° f dl Kp — > — ° / dl K(p, 

where <p is the Hilbert transform of ip on the boundary of the domain D, 
with the property di<p\ dD = —d n ip\ dL) . The normal derivative is taken in the 
direction inside D. Once cp is harmonic p coincides with the definition pip . 
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The electric-magnetic transformation makes (p a pseudo-scalar. As a result 
the action reads 

A = -!- [ d 2 x CVip) 2 - f dl K<p, dilute phase, a < 0. (20) 

J D 2ir J dD 

We note that in the dilute phase the action is real and its saddle point 
lies in the real fields. The connection of the field ip to the critical curves in 
the dilute phase will be discussed below. 

For simplicity, we start with the theory in the upper half plane. The 
components of the stress-energy tensor now read (cf. (j!2p ) 

T(z) = ~(dcj)) 2 + iV2a d 2 <j), f(z) = -^(54>) 2 -iV2a o d 2 0. (21) 

The condition T(x) = T(x) now glues the holomorphic and antiholomorphic 
sectors on the real axis as (f>(x) = —<j)(x) which means (p(z, z) = (j)(z) — 4>(z), 
<p = —icj)(z) — i(p(z) and therefore the boundary condition on the real axis is 
Dirichlet: 

<M,= = °- ( 22 ) 

Unlike in the dense phase, (ft is complex conjugate of 4>. 

The holomorphic and antiholomorphic weig hts of 0< e ' m ) in the dilute 
phase are (cf. (fTTj) ) 

h(e, m) = (e + m)(e + m — 2ao), h(e, m) = (e — m)(e — m + 2ao) 

or 

h = a(a — 2ao), h = 6t{a + 2ao) 

in the notation (|15p . Since in the dilute phase «o < 0, the charge which 
corresponds to a weight h and vanishes with h is 

«h = «o + \j oiQ + h, k ^ 4 (23) 

A spinless vertex operator in the dilute phase is either purely magnetic e = 0, 
or it has a magnetic charge m = oiq and an arbitrary electric charge. 

The stress-energy tensor f|21 1) shows that under a conformal transforma- 
tion z — ► f(z) 

<p(z,z) -> (p(f(z),f(z)) - 2v / 2a arg/'(z), 
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3D 

Figure 7: Boundary condition ([25]) . i? = — arg f' D | QD . 



which preserves the reality of (p. We observe that the flow lines of the vector 
field 



e 2V2 ao * (24) 

are conformally invariant since (|24p has conformal spin h—h = — 1. li (p were 
a smooth function the flow lines would repulse. This suggests identifying 
them with the critical curves in the dilute phase. The identification of flow 
lines (|24p for the Gaussian field with Dirichlet boundary conditions with 
conformally invariant curves in SLE was done in [39]. 

In the dilute phase the boundary of the domain is a flow line. The 
boundary condition (I22p transforms with the transformation of the domain. 
If /d(-z) maps the domain D onto the upper half plane (a domain with a 
straight boundary) the boundary condition in D is (Fig. [7]) 

v\dD = - 2 ^ a oargf' D \ dD . (25) 

At one point on the boundary ip undergoes a jump 4\/2ao 7r ) which means 
that a magnetic charge 2ao is inserted there. In the upper half plane this 
magnetic charge is pushed out to infinity. 
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In the upper half plane the Dirichlet boundary condition allows us to 
write the vertex operator (fT4"|) as (cf. eq. (fl6j) ) 

( ' e ' m \z,z) = e *v / 2(e+m)0(2:) e iv / 2(m-e)0(z)^ 

which shows that the image of an electric charge changes sign and the image 
of the magnetic charge does not. 

The conformal weights and holomorphic charges of the curve-creating 
operators are found by the continuation of the corresponding formulas from 
the dense phase. That is, the bulk operators giving rise to n critical dilute 
curves are ^0^/2 an d V'l.n+i, but with k < 4 in the definitions (p~9|) of their 
weights. We summarize the operators which create n curves in the bulk and 
on the boundary in both phases in the following table. 





n curves in bulk 


n curves on boundary 


r, s 


e, m 


explicit 


r, s 


a 


explicit 


Dense 


- 

u, 2 






l,n+l 




-i-2=a_<A 
e V2 v 


Dilute 


- 

u, 2 


=Ff a_,a 




l,ra + l 




-i-^Q_0 



4.3 Differential equations for the boundary curve-creating 
operator 

Due to the conformal invariance any correlation function which includes a 
boundary curve-creating operator satisfies a differential equation which we 
derive below. The correlation functions then can be found as solutions of 
this equation. This applies to both the dense and the dilute phases. 

The stress-energy tensor (|12|) or (|21|) is a field which generates trans- 
formations of primary fields under changes of geometry. The change of the 
upper half plane into that with a puncture at wq is produced by an infinites- 
imal transformation conformal everywhere except Wq: 

f( z ) = z + e (—!— + —L-), e^O. 

V z — Wo z — Wo-' 

This transformation affects correlation functions of primary fields in the 
usual way via the primary transformation law. A correlation function in the 
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new geometry is 



k k / \ t k 

M\w 



i=l i=l i=l 

i + « E (tt^-w - + {»» - *>)] ( n v« i . 



-A - io ) 2 Zi-^o 



i=l 



The same change of the domain is by definition produced by the insertion 
of the stress-energy tensor T(wo) + T(wq) into the correlation function (the 
conformal Ward identity). 

We now fix wq to be a real number x. On the real axis T(x) = T{x) and 
the Ward identity reads 

(T(x)0(zi,. . . z k ))u = C- 2 (0(zi, ■ ■ ■ z k ))m, 
where 0(z x , . . . z k ) = V 011 ^) . . . V ak {z k ) and 

^-E(^-^^)- <*> 

1=1 ' 

We remark that due to translational invariance the action of d x + Yli=i 
on the correlation function is zero so that C-i = d x . 

The product of fields inside a correlation function is not normal-ordered. 
We can also insert the stress-energy tensor at the same point as one of the 
fields using normal ordering: 

(:T(x)V a (x) : 0(zi, . . . z k )) m = C- 2 {V a (x)0{ Zl , . . . z k )) m (27) 

with the same definition of £-2- For special values of a this leads to a 
differential equation [2j [3] . 

Using the series expansion V a (x) = Yl^Lo ^^"^ -4> n (x) ■ and the Wick 
theorem we find the normal-ordered product 

:d<p(x)d(f)(x)V a (x): = :d4>(x)d(j){x)e iV2a ^ : -i2y/2a :d 2 cp{x)e iV2a ^ :, 
:d 2 (^{x)V a {x): =:d 2 ^(x)e iV2a ^ x) : . 

The second term in the first line comes from only one of d<f){x) being con- 
tracted with the exponential. We then find 

:T{x)V a {x) :=: (-^d(j)(x)d^(x) + iV2(a + a)d 2 (f>{x))e iV2a ^ x) : . 



22 



On the other hand, a simple differentiation gives 

d 2 V a {x) =:{-2a 2 d(l){x)d(t>{x) + iV2ad 2 0{x))e iV2a</,{x) : . 

These two expressions are the same up to a prefactor provided 4a(a + qo) = 
1. The two solutions of this equation are ai j2 = \ j \[k and 02,1 = — vW^- 
They determine operators degenerate on level 2: = e 1 ^ 2 " 1 - 2 ^ and tp2,i = 
e iV2a 2 ,i4>. We then write 

: (^8 2 - T(x))Mx) : = 0, : (V - T(x))^,i(x) : = 0. 

These are operator equations, meaning that they can be inserted into any 
correlation function leading, for example, to 

(: {^d 2 x - r(x))^i, 2 (x) : 0(3i, • • • z fc )>H = 0. 

If all fields are primary we combine this result with the Ward identity (|27p 
and obtain the following differential equations: 

(^C 2 _t - C-2)(Mx)0( Zl , . . . z k )) m = 0, (28) 

{-Cli - £-2)^2,1(^)0(21, ■ • • z k ))n = 0. 

The first of these equations is relevant for calculations of correlation func- 
tions containing the boundary curve-creating operator ^1,2- 

5 Stochastic (Schramm-) Loewner evolution (SLE) 

In the previous sections correlation functions of boundary operators were 
interpreted as averaging over the shapes of critical curves. The shape of 
a curve fluctuates with some probability measure i.e. there is a certain 
statistical weight associated with each possible shape. The SLE approach 
[8] determines this measure. If a curve is conditioned to a certain shape 
objects of field theory explicitly depend on it. 

In this section we briefly review SLE ([8]-[19]) and present the Bose field, 
its current, the vertex operators and the stress-energy tensor explicitly in 
terms of the shape of the curve. We restrict the discussion to the dilute case 
k € 4. 
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5.1 Loewner equation 

Consider the upper half-plane EI in which an arbitrary non-self-intersecting 
curve 7 starting from the origin is drawn. Let the curve be parameterized by 
t ^ called "time" . This setting defines a continuous family of slit domains 
each member of the family being the upper half plane with a part of 7 up 
to some value t removed: Ht = H\7t. Each slit domain Mt is topologically 
equivalent to the upper half plane since 7 does not have self-intersections. 
We define the time-dependent conformal map wt(z) of the slit domain to 
the upper half plane: — > H. To make the definition unique we specify 
the behavior of wt(z) at infinity: 

w t (z) = z- & + — + ..., z^oo, (29) 

where £t is a real function of time uniquely determined by the curve jf The 
time-evolution of the map is given by the Loewner equation [40] 

2 

d t w t (z) = — wt(2) L_ = 2, &L_o = °- (30) 

At each z this equation is valid up to the time t when wt(z) = 0. When 
this happens the point z no longer belongs to the slit domain (it has been 
swallowed by the curve). For a simple curve this happens when z is hit by 
the curve: jt = z. Hence, the Loewner function wt(z) maps the tip of the 
growing curve onto the origin. 

It is obvious that in order for 7 to be a continuous curve without branch- 
ing the forcing ^ should be a continuous function. Actually, certain types of 
forcing produce curves with double points: 7 can self-touch. On the event 
of self-touching a portion of the plane that has become encircled by 7, is 
removed from M together with the [0, t] portion of the curve, thus ensur- 
ing that the topology of does not change. The applicability range of 
the Loewner equation, therefore, extends beyond the originally planned de- 
scription of one-dimensional objects, e.g. curves. In general, the objects that 
grow according to the Loewner equation (the unions of swallowed points) 
are called hulls. 

The Loewner equation (|30|) with a specified forcing defines a curve or 
a hull in the upper half plane growing with time, the process known as 
the chordal Loewner evolution. Its basic property is that the final point 
of the evolution is lim^oo 74 = 00 i.e. a point on the boundary of the 
domain H. One can define the chordal Loewner evolution in any domain D 
topologically equivalent to H by a simple composition of maps. In this case 
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by the Riemann theorem there exists a conformal mapping f(z) : D — > H, 
and wt(f(z)) defines a growing curve 74 in the domain D: it maps D\jt — > H. 
The final point of this evolution / _1 (oo) necessarily lies on the boundary of 
D. 

5.2 Stochastic Loewner equation 

The Loewner evolution can be used to describe any curve which starts from 
the origin and goes to infinity in the upper half-plane. We would like to 
apply it to a conformally invariant curve produced in CFT by the insertion 
of operators tpip a t and at 00. This dictates the necessary properties of 
the forcing 

Since 7 in a statistical system is a random curve the Loewner forcing ^ 
should be a stochastic process. This turns the Loewner equation (|30f) into 
a stochastic differential equation |41| W2\ which has the form of a Langevin 
equation. In the rigorous sense it should be written in terms of stochastic 
differentials because £t may (and will) be nowhere-differentiable: 

dw t (z) = w t +dt(z) ~ w t {z) = 2< f - d£ t , d£ t = €t+dt - €t- (31) 

w t {z) 

The r.h.s. is evaluated at the time t (the Ito convention for stochastic anal- 
ysis |41j). This remark is important when we deal with Markov processes: 
their behavior at any time is completely independent of their earlier behav- 
ior. In the following we will write stochastic differential equations with usual 
time-derivatives, as in (|30p . always understanding them in the ltd sense. 

The properties that specify the driving function are obtained from 
the following simple considerations. First of all, a critical curve 7 in a 
translationally-invariant system has no special preference of going to the 
left or to the right. The £t should have a distribution that is symmetric 
around zero leading to the requirement 

6 = -£t in law. (32) 

By this we mean the statistical identity: random quantities on both sides of 
the equation are taken from the same distribution but, in general, represent 
different realizations. 

Next, the averaging over the ensemble of shapes of 7 can be done in 
steps: first we can fix the shape of the curve up to a time t i.e. 7[0,t], and 
average over the remaining part j[t, 00]. In this averaging 7[0, t] is seen as 
a part of the domain boundary whose shape will be averaged later. With 
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Figure 8: The image of the [t, t'] part of the curve on the left created by the 
mapping wt{z) is statistically identical to the [0,t"] part of the same curve 
on the right. 



a fixed 7[0,i] the remaining part j[t, oo] is a conformally invariant curve 
in the domain H^. By the conformal invariance we can map Mt onto HI by 
wt(z) and the image of j[t, oo] should be statistically identical to the original 
whole curve 7(0,00] (see Fig. [8]). In particular, the image of any segment 
7[t,t'] (here t' > t) should be statistically identical to a segment 7(0, t"] of 
the original curve: 

7[0, t"] = w t (j[t, t']) in law. 

The fact that we can perform the averaging in steps means that the r.h.s. 
and the l.h.s. are statistically independent. The application of Wf maps the 
curve in the l.h.s. into the real axis, and the same is done with the curve 
7[i, t'] by the application of w#. Therefore, w t '(z) = Wf (wt(z)) in law. At 
large z we can replace the Loewner maps with their asymptotic form (|29p 
which allows us to deduce t" = t' — t and arrive to two equivalent statements 

w t '(z) = w t '-t(w t (z)) inlaw, (33) 
it' =it + it'-t in law, (34) 

where the r.h.s. and the l.h.s. are statistically independent in both formulas. 
Instead of eq. ([33]) we may equally well write 

w t i(z) = w t {w t >-t{z)) inlaw. (35) 

The only continuous stochastic process & with stationary independent 
increments as in eq. (|34|) and the symmetry (|32|) is a driftless Brownian mo- 
tion. Thus, it follows from the conformal invariance and continuity of the 
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curve that it should be described by stochastic Loewner evolution (SLE). It 
is defined by eq. (|30p where the forcing ^ is a Brownian motion. The re- 
quirement of continuity is important here: it ensures that the curve does not 
branch and can be chronologically ordered in a unique way. Hence, the prob- 
ability measure on a conformally invariant curve starting from the boundary 
is related to the Brownian motion via the Loewner equation. There is only 
one free parameter here: the strength of the noise k defined by 

where we use the notation -< . . . >- for averages over realizations of the 
Brownian motion 

Different values of k must give rise to conformally invariant curves in dif- 
ferent statistical systems. Actually, the curve 7 is simple i.e. never touches 
the boundary of the domain (or itself, by conformal invariance), if and only 
if k ^ 4 [31] as in Fig. [H This is consistent with parametrization of systems 
with dilute critical curves by k < 4 in Section [H At k > 4 the Loewner 
transformation wt(z) maps the external perimeter of 74 onto the real axis. 

Next we turn to the formula (|35p which we rewrite using the sign o for 
composition of maps as w^ 1 o Wf = Wtf—t in law. Note that w^ 1 o Wf maps 
Mf onto Mf. Consider now an infinitesimal time increment: dt = if — t, for 
which we have w^f 1 o wt+dt = Wdt m law. The infinitesimal Loewner map in 
the r.h.s. is given by (see eq. (|31|)) 

Wdt (z) = z + — -di t . (36) 

z 

Then the statistically equivalent infinitesimal map w^ 1 o Wt+dt from Mt+dt 
to M t is given by the same equation: 

Wt 1 (w t+ dt(z)) = z + — - d£ t - (37) 

We stress, however, that the equivalence of these maps is only statistical 
and the Brownian motions £j in the two formulas (|36l |37|) are two different 
realizations drawn from the same statistical ensemble. Treating them as the 
same function would lead to an inconsistency because no function £t (apart 
from a linear function at) satisfies the requirement ()34|) in the exact, not 
statistical, sense. 

5.3 Martingales and correlation functions 

If a stochastic process 0(t) is a martingale, its average does not depend on 
time. Below we explain how a martingale 0(t) generated in SLE corresponds 
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to a CFT holomorphic operator O such that [lOj 



(^l,2(0)Vl,2(oo))H 

The identification is up to a constant Co which depends on aQ. We do not 
discuss it here. 

Comparison with the two-step averaging ([2]) in which, as we now know, 
the curve-creating operator is B = tfti 2, suggests the identification 



0(t) ~ (0* (00)) 



it ■ 



We will demonstrate the correspondence in the case of primary fields. 

We will further explore the relation between CFT and SLE and show 
that at t — > 00 the SLE fields 0(t) satisfy the operator algebra of their 
counterparts in CFT. This is to be expected. The size of 74 introduces a 
length scale ~ yi. When it diverges at t — > 00, the local fields in SLE 
fluctuate in a scale-invariant way. In the following discussion we consider 
only the dilute phase k ^ 4. 



5.3.1 Primary fields in SLE 

The Loewner equation can be considered in any domain. The fields related 
to the evolving curve transform under the conformal map / which maps one 
domain onto another. We call an SLE primary field a field which evolves 
under the Loewner equation and transforms as Oh{z) — > f' h (z)Oh{f{z)) 
under a change of the domain. In particular, we can consider a pair of 
domains H and Eh, related by conformal map Wt{z) in which case 

O h (z,t) = w' t (z) h O h (w t (z),0). (39) 

Alternatively, we can use the pair of domains Hh, and Wt+dt related by 
an infinitesimal conformal transformation w^ 1 (wt+dt( z )) in which case 

O h (z,t + dt) = (d z Wt l (w t+dt (z))) h O h (wt l (w t+dt (z)),t). 

Now we use the stochastic nature of SLE: the statistical identity of w^ 1 o 
wt+dt and Wdt and the explicit expression (f37|) enables us to write 



O h (z,t + dt) = w' dt (z) h O h (w d t(z),t) = 
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Expanding the r.h.s. to the first order in dt and to the second order in d£t 
produces the evolution equation for Oh(z,t): 

d t O h (z,t) = (^Cl 1 -2C- 2 + iC- 1 )o h (z,t), (40) 

where we used the notation (|26p with k = 1 and the position of the curve- 
creating operator x set to 0. The derivation is extended to a product of k 
primary fields, when operators £_2 act on k variables. 

The two ways of representing the stochastic evolution of Oh(z,t) in the 
time-derivative of eq. ([39]) and in eq. (f4T)|) are equivalent but cannot be 
obtained from each other by a simple substitution since they contain different 
realizations of Brownian motion. 

Since we understand stochastic equations in the Ito sense, the differential 
of the Brownian motion ^t+dt ~ £t in the last term of eq. (|40p is multiplied 
by d z Oh(z,t) evaluated at time t. Therefore, by the Markov property of 
the Brownian motion this term vanishes upon taking the average. The 
expectation value of a primary field then satisfies the Feynman-Kac partial 
differential equation 

d t <O h (z, ty = (^Cl 1 - 2£_ 2 ) -<O h (z, t)y , (41) 

and similarly for a product of k fields. The evolution equation for a primary 
field is obtained solely from its transformation law under an infinitesimal 
conformal map. The same method allows to obtain the evolution equation 
for any field with a specified transformation law, not necessarily primary. 
An example of a primary field is 

O h (z,t)=t a ^w' t (z) h , (42) 

using the notation (fT9"j) and (f2"3"j) . As follows from the scaling law of SLE: 

wAz) = -w a 2 t (za), in law, 
a 

the solution of the equation (|4ip for a real z = x in this case must have the 
form t ai ' 2ah F(x/\/t) with the finite asymptote x 2ai ' 20th at t — > oo. In this 
limit the fluctuations of Oh become scale-invariant and its average corre- 
sponds to a correlation function in CFT. There local fields have a singular 
short distance expansion. We now comment on how it arises in the SLE 
approach. 

Consider the product of two boundary fields t ai - 2 ^ ah i +ah ^w' t (xi) hl w' t (x2) h2 ■ 
Taking x\ = X2 we obtain t ai < 2 ( ah i +ah '^w' t (xi) hl+h2 , which diverges at t — > 
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oo due to convexity of ah- The properly normalized field t ai,2a( ~ h ^ +h ^w' t {xi) hl+h2 
is different. This is an example of the short distance singularity It closely 
parallels the origin of such singularities in field theory. 

Since dt <Oh{z,t)>- vanishes at t — > oo, the Feynman-Kac equation 
(|4ip in this limit becomes identical to eq. ([58]) . Therefore, the SLE average 
-<Oh(z, oo)>- and the CFT correlation function (ipi,2(0)V ah (z)^(oo))m obey 
the same equation. However, this fact alone does not allow to identify them 
as the following simple remark shows. The CFT correlation function is a 
holomorphic function z 2ai - 2ah with a branch point at the origin. It acquires 
a phase in going from a real z to —z. However, the SLE field t ai > 2ah w' t (z) h 
does not since w' t is real on the real axis. 

The proper choice of the primary field of weight h ia^| 

V° h (z, t) = w t {z) 2ai < 2ah w' t {z) h . (43) 

Its expectation value is holomorphic and can be identified with the holomor- 
phic CFT correlation function (ip h2 (0)V ah (z)^(oo)) = z 2ai ' 2ah . Indeed, its 
initial value is V ah (z, t)| t _ = z 2ai ' 2ah and as we will show momentarily, its 
average does not depend on time. 

Using the Loewner equation and the Ito formula for a product of two 
stochastic processes X t and Y t : 

d(X t Y t ) = X t dY t + Y t dXt + dXtdYt (44) 

with X t = w t (z) 2ai - 2ah and Y t = w' t {z) h we find 

d t v a "(z,t) = -^f^v a HzM, 

w t {z) 

which yields d t -<V ah (z,t)y = i.e. V ah (z,t) = w t (z) 2ai < 2ah w' t (z) h is 
a martingale. We note that a\p, = and thus a\^t is a Brownian 

motion Bt with variance 1. 

We conclude that -<V ah (z, t)y satisfies the stationary equation which is 
identical to (1251): 



-CU-C^^V ah (z,t)y = Q. 

at any time and not only at t — > oo. The same equation holds for a string 
of k primary fields. 

In the remainder of this section we will further illustrate the correspon- 
dence between SLE and Bose field. 



2 This formula is understood up to a multiplicative constant related to the constant 
C o {a ) in 438}. 
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5.3.2 Bose field in SLE 



We now will discuss the SLE Bose field 4>(z,t), a martingale which corre- 
sponds to the holomorphic part of the Bose field discussed in Sec. 14,21 

4>{z,t)~ (v>i, 2 (o)<K^(oo)> 7t . 

It follows that 4>(z, t) has the transformation law (|13|) and is a martingale. 
This fixes it to be [15] 

iTZ 

4>{z,t) = iV2a log w' t (z) - — log w t {z), TZ = 

up to an additive constant (cf. footnote [2]). Indeed, it follows from the 
Loewner equation that dcf>(z, t) = d£,t and therefore the average of 4>(z, t) 

equals its initial value -<4>(z, t)y = — ^ log z. The latter is just the holomor- 
phic part of the solution of the Laplace equation with a step-like boundary 
condition ip(x < 0) = ttTZ, ip(x > 0) = 0. 

The real field cp(z, z, t) = (p(z, t) + cj)(z, t) is compactified with radius TZ. 
If we choose a point z 7 to lie on the slit 7 the value of 99 is given by the 
direction $ of the tangential vectoJl: 

(pL(zy) = 2V2a Q ti + ttTZ, v? fl (z 7 ) = 2^00(1? - vr), (45) 

where the additive constants correspond to the two sides of the slit (Fig. [9]). 
The Bose field then measures the winding angle of 7. The tangential vector 
of 7 is e i ¥'fl/ 2 v / 2ao 5 as i n (Jjzjj. The identification of the Bose field with an 
angle emphasizes that in the dilute phase it is a pseudo-scalar. The field 
fluctuates with 7 but its jump across 7 is fixed: 

(<PL ~ VR% = <K + 2V2«o) = yi^r. 

On the real axis ip obeys the Dirichlet boundary condition: <p(x, t) = for 
x > and <p(x, t) = ttTZ for x < 0. 

The case n = 4 (i.e. «o = 0) is special: the values of the field on jt then 
are fixed to be and ttTZ as on the real axis. 

The evolution equation for the Bose field follows from its transformation 
law and the Loewner equation similarly to the previous section. It is given 
by 

d t <P(z,t) = g/^ - 2£_ 2 + " (46) 

3 In fact, 7 is a fractal curve. It is understood that 1? is measured not on 7 but nearby. 
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Figure 9: The SLE field (p on 7 measures the winding angle •& (see eq. (|45j) ). 
The tangential component of the current dup on 7 measures the curvature. 



with h = in the definitions fl26f> of £'s. Since (p(z, t) is a martingale its 
average -<4>(z, t)y = — ^ log z solves the equation 

(^_ 1 -£_ 2 )^(,,^ -^ = 0. (47) 

This is the same equation as (ipi t 2(0)4>(z)^(oo)) satisfies. 

We now comment on how the operator algebra of the Bose field arises in 
the SLE approach. With the help of the Ito formula (|44p we find from the 
Loewner equation the evolution of the product of two Bose fields: 

d(4>(z 1 ,t)(p(z 2 ,t)) = d(p(zi,t)(p(z2,t) + 4>(zi,t)d(j)(z2,t) + d(p(zi,t)d(j)(z2,t). 

One can check that the last term (drift) is 

d<j)(z 1 ,t)d(j)(z2,t) = t~V~1 — \ = d\og(w t (zi) - w t (z 2 )), 

wt{zi)w t (z 2 ) 

and therefore 4>(zi,t)(j)(z2,t) is not a martingale. The martingale reads 

<f>(z u t)<f>(z2,t) - log(iwt(zi) - w t {z 2 )). (48) 

This process is identified with the product of operators in CFT. The sub- 
tracted term here is non-local. However, at t — > 00 it vanishes — a manifes- 
tation of the fact that only in this limit local fields of SLE fluctuate as local 
fields of CFT. 
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In the limit z\ — > zi the subtracted term becomes \og{wt(z\) — Wt{z2)) = 
log w' t {z{) + log(zi — z%). We define the normal-ordered product by subtract- 
ing the logarithmic divergence from (|48D : 



:<t> 2 (z,t):= <j> 2 (z,t) - log^). 

The ltd formula (|44p applied to a product of more than two bosons be- 
comes the Wick theorem. The normal-ordered exponential is found up to a 
multiplicative constant as 

which is just the definition (|4"3"j) of i). 
The holomorphic current of the Bose field is 

% r- w'J(z) IZw'Az) 

j(z,t) = 84>{z,t) = iV2a ^rri - *T-TT 

w' t (z) 2 io t (z) 

Since cp is a pseudo-scalar, the current is an axial vector. On the curve 7 its 
tangential component is proportional to the geodesic curvature of 7: 

ji\ y = 2V2a K T 

This is consistent with the action (|2U|) . At the tip of the curve 74 the current 
has a vortex with a charge proportional to «o- 

The evolution of j(z,t) and the equation for -<j(z, t) y are found by 
differentiating eqq. (|46l47p . Differentiating (j4"5|) and taking zi — ► Z2 we get 
the normal ordered square of the current 

■f{z,t):= j 2 {z,t) - -{w t ,z}, 


where {w,z} is the Schwarz derivative: 

, , w"'(z) 3fw"(z)^2 
{w,zj - 



w'{z) 2\w'(z) 

Finally, we can define the holomorphic stress-energy tensor in SLE as a 
martingale with the proper transformation law: 

z - f(z) : T(z,t) - f'(z) 2 T(f(z),t) + ^{f,z}, 
with c = 1 — 24«q. This fixes it to be 

T{z,t) = -- :d<p(z,t)d<p(z,t): +iV2a o d 2 0(z,t), 
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which explicitly reads 

r^ f % c , -, ( w'Jz)\ 2 , 6 — K 

T(z,t) = -{w t ,z} + h li2 [^y-) , h l , 2 = ——. 
12 \wt(z) ) 2k 

The combination T(z,t)dz 2 can be used to determine the probability that 
7 intersects a segment of size dz centered at z \43\ [4"4"j . 

The averages of the current and the stress-energy tensor are 

i\J2a\ 2 fei o 

z z z 
We note that the current and the stress-energy tensor measure the charge 
and the weight of the curve-creating operator ipi^- 

We thus have seen that the stochastic Loewner evolution allows to write 
the objects of field theory explicitly in terms of the shape of the curve. These 
objects can be used to find distributions of various characteristics of curves. 
As an application we compute in the next section the harmonic measure of 
critical curves for k ^ 4. 



6 Harmonic measure of critical curves 

In this section we apply the results of the previous sections to compute the 
multifractal spectrum of the harmonic measure of dilute (k ^ 4) critical 
curves. In the dense phase k > 4 these results apply to external perimeters 
of critical curves. 

Harmonic measure quantifies geometry of complicated plane domains 
[45], [46] . We start with basic definitions. 

6.1 Harmonic measure of critical curves 

Consider a closed simple curve 7. One can imagine that it is made of a 
conducting material and carries a total electric charge one. The harmonic 
measure of any part of 7 is defined as the electric charge of this part. In 
what follows we will pick a point of interest zq on the curve and consider a 
disc of a small radius r <C R centered at zq. It surrounds a small part of 

7 and we define (i(zo,r) to be the harmonic measure of this part i.e. the 
electric charge in it. 

Then we consider the moments 

N 

M h = J2^r) h , (49) 
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Figure 10: A curve covered by discs. 



where h is a real power and N is the number of discs needed to cover 7 as in 
Fig. [TUl As the radius r gets smaller and the number of discs A?" gets larger, 
these moments scale as 



where L is the typical size of the curve. 

The function r(h) is called the multifractal spectrum of the curve 7. 
This function encodes a lot of information of the curve 7. It also has some 
simple properties. First of all, since all < fi(zi,r) ^ 1, the moments 
Mh are well defined for any real h and the function r(h) is non-decreasing: 
r(h) ^ T~(h') for any h < h! . Secondly, if h = 1 the sum in (|49p is equal to 
the total charge of the cluster and therefore does not scale with r, producing 
the normalization condition r(l) = 0. Thirdly, if we set h = 0, Mq is simply 
the number of discs of radius r necessary to cover the curve 7 so that by 
definition the fractal (Hausdorff) dimension of 7 is df = — r(0). 

If the curve 7 were smooth we would have a simple relation r{h) = h — 1. 
For a fractal curve one defines the anomalous exponents 5(h) by 



Also, the generalized multifractal dimensions of a fractal curve 7 are defined 



as D(h) = r(h)/(h - 1) (so that D(0) = d f ). A theorem [IS [48] states that 



If the curve 7 is random, so are its moments , and we can average them 
over the ensemble of curves. It is natural to assume that the summation 




r(7i) =h-l + d(h). 



D(1)=t'(1) = 1. 
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over the points Zi in eq. (j4"9j) is equivalent to the ensemble average. Hence 
we can write 

<M h y = iV^u(r)V ~ (^) T(0) ^(r)V , 

where now the harmonic measure fj,(r) is measured at any point z on 7 and 
-<...>- indicates the ensemble averaging. We define the local multifractal 
exponent r(h) = r(h) — r(0) at a point z by 

^( P )V ~ ( I ) . (50) 

The ensemble of critical curves is suitable for multifractal analysis. There 
are a few generalizations of the simple closed curve considered above. First of 
all, the curve 7 need not be closed or stay away from the system boundaries. 
If 7 touches a boundary we can supplement it with its mirror image across 
this boundary. The definition of fi(z, r) can be naturally extended to the 
cases when z is the endpoint of n critical curves on the boundary or in the 
bulk. If n is even, the latter case can be also seen as n/2 critical curves 
passing through z. In particular, n = 2 corresponds to the situation of a 
single curve in the bulk considered above. 

When z is the endpoint of n critical curves on the boundary or in the 
bulk we define the corresponding scaling exponents similarly to eq. f|50j) : 

M^rfy ~r h+ ^ h \ M^rfy ~r ft <("). (51) 

In the case of a single curve we will drop the superscript, for example, 
A(h) = AW(/i). 

These exponents were first obtained by means of quantum gravity in 
[U [5] . For a critical system parameterized by k the results read 

k - 4 + y/(n - 4) 2 + 16Kb _ VI - c + 24/t - VT^c 
2~k V25^- ' 

nA(h), (53) 

~\ + ( T6 + ^) ( K " 4 + " 4) 2 + 16«/0 . (54) 

Remarkably, A(/i) is the gravitationally dressed dimension h as given by 
the KPZ formula of 2D quantum gravity [7j. Connections of this formula to 
SLE were discussed in [35] . Starting in the next section we will show how to 
obtain these exponents formulating the c ^ 1 CFT via the Bose field, where 
they are also written in a more transparent way. 



A(h) 
A( n \h) 
A buikC0 
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Figure 11: The uniformizing conformal maps for various cases considered. 
The dots denote that points where the electric field is measured. 

6.2 Harmonic measure and uniformizing conformal maps 

Let us consider a conformal map w(z) of the exterior of 7 to a standard 
domain. Usually we will choose the upper half plane but sometimes the 
exterior of a unit circle is more convenient. We choose the map so that 
the point of interest z on the curve is mapped onto itself, choose it to be 
the origin z = and normalize the conformal map so that w(oo) = 00 and 
w'(oo) = const. In the case of the upper half plane we choose w'(oo) = 1 as 
we did in SLE (see Section [5]), while for the unit circle w'(oo) = p" 1 where 
p is the conformal radius of 7. 

The scaling of w(z) near the origin is directly related to that of the 
harmonic measure. Indeed, since fi(0, r) is the charge inside the disc of radius 
r by the Gauss theorem it is equal to the flux of the electric field through 
the boundary of this disc. This, in turn, should scale as the circumference 
of the disc times a typical absolute value of the electric field at a distance r 
from the origin i.e. |u/(r)|. This leads to a scaling relation 

n(0, r) ~ r\w'(r)\, 

which allows to rewrite the definitions (I5ip as 

-V(r)|V ~r A(n) W, -V(r)|V ~r A b">). 

The relation between the scaling of the harmonic measure and the deriva- 
tive of a uniformizing map allows further generalizations. Namely, we can 
measure the electric field in more than one point. Consider n critical curves 
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emanating from the origin. Close to the origin the curves divide the plane 
into n sectors if the origin lies in the bulk and n + 1 sectors if it lies on the 
boundary. Then we can study objects like 

■<H\w\ Zi )\ h ^, 

i 

where no two z^s lie in the same sector. The case when the electric field is 
not measured in some sectors is done by setting hi = in them. We will 
see how to express these quantities as CFT correlation functions. In the 
case when Zj's are all at the distance r from the origin these averages scale 
as r A(n) Oi.---'>.7i+i) anc j r A Lik( h i'"' h ™) with the higher multifractal exponents 
related to the multifractal exponents ([521 - as [5] 

n+l n+1 

A^(ht, . . . h n+1 ) = A(n) (hi) + \ XJ A (k) A (M, (55) 

i=l i<j 
n n 

A$L(hu- ■■*») = £ A bl(^) + \ E A^A^). (56) 

i=l i<j 

In the case when z^s are independent variables the answer is more compli- 
cated. However, in the case of a single curve at the boundary it reads: 

^\w\x 1 )\ h '\w'(x 2 )\ h ^ ~ H^H^IX! - X2 |A(fci,h 3 )-A(h 1 )-A(h 3 ) > 

(57) 



6.3 Calculation of multifractal spectrum 

Here we obtain formulas (f52l - l54l l55l l56j) . We will consider in detail the 
simplest case of harmonic measure of a single curve near the boundary. The 
other cases contain a few technical differences. The calculation is based on 
the two-step averaging of Sec. El 



6.3.1 Boundary multifractal exponent from CFT 

All averages over the ensemble of critical curves are some correlation func- 
tions of the Bose field. Averages of w' h are related to correlation functions of 
primary fields Oh in the presence of curve-creating operators. We start with 
the case of a single curve starting from the boundary, when this operator is 
V>i,2 and find the boundary multifractal exponent A(h). 
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We use the two-step averaging result ([2]) where we write B = ipi^ and 
choose O to be a boundary vertex operator Oh(r) = e l ^ ah ^ r h We uni- 
formize the exterior of 7 by mapping it to the upper half plane. Being a 
primary operator of weight h, Oh(r) transforms as Oh — ► w'(r) h Oh{w(r)) 
while ^(00) does not transform because of the normalization of w(z) at in- 
finity. The transformation relates the correlation function in the exterior of 
7 to a correlation function in the upper half plane: 

(0 & (r)tt(oo)> 7 = w'(r) h (O h (w(r))t>(oo)) m . (58) 

Note that the correlation function in the r.h.s. has no singular dependence 
on r, hence, as long as we stay in the limit r <C \L\ it is a constant prefactor. 
We then obtain a scaling relational 

-fti/(r)V ~ (0 h (r)Vi,2(O)tf(oo))H. 

The primary field ^(oo) should be chosen such as to satisfy the physical 
condition A(0) = 0. This formula establishes the relation between the mul- 
tifractal spectrum of critical curves and the correlation functions. Explicitly, 
the correlation function reads ^V2a h 4>(r) e iV2a h2 m ^(^-^-^1,2)^00)^ 

where 0, r, 00 are three points on the real axis. It remains to find how it 
scales with r. 

To find the r-dependence of this correlation function, we use the fusion 
of vertex operators (I18j) and obtain the result (152p written in a compact and 
suggestive form: 

A(h) = 2ai i2 ah, ah = "o + yJo$ + h. (59) 

An immediate generalization to the statistics of the harmonic measure 
of n curves reaching the system boundary at the same point is obtained by 
replacing -01,2 — ► ipi,n+i- Since ai, n +i = —na-/2 = na\ t 2, this immediately 
leads to 

A^(h)=2a hn+l a h = nA(h), 

which is the same as eq. ([531) . 

We also remark that nothing compels us to measure the electric field 
on the real axis. Instead of the boundary field Oh{r) we could take a bulk 
primary field Oy{z, z) where \z\ = r. The weight h! should be chosen such 
that when the holomorphic part Oy(z) is fused with its image Oh'(z) the 
boundary field Oh(^y^) with weight h is recovered. 

4 We drop the absolute value since on the real axis w' is positive. 
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6.3.2 Boundary multifractal exponent from SLE 

For comparison, we compute the multifractal spectrum of a single curve in 
the SLE approach (see e.g. [Tl"]). 

In SLE we must compute ^w' t (r) h >- in the limit t — > oo since the size 
of the curve L ~ \ft. With the choice of the primary field as in (|42p the eq. 
(|4"Tj) at t —> oo reads 



The solution is a power law -<w' t (r) h >- ~ r 2 " 1 - 20 ' 1 , in agreement with (|59l) . 

6.3.3 Calculation of higher boundary multifractal exponents 

We consider n non-intersecting critical curves growing from the origin on the 
boundary (the real axis). It will be convenient to assume that the curves end 
somewhere in the bulk thus forming a boundary star (e.g. the third picture 
in Fig. [TT]) . Let w(z) be the conformal map of the slit domain formed by 
the star to the upper half plane normalized at infinity. 
We want to find the scaling of the average 



where Zj's are all close to the origin, no two of them lying in the same 
sector. The latter condition will be automatically satisfied in the subsequent 
calculation due to the following: if in a particular realization some two points 
Zi and Zj happen to be in the same sector, then w{zi) — w(zj) — > as the 
size of the star L — > oo. 

Since n curves starting from the origin on the boundary are produced 
by the operator ipi tU+ i(0) we now consider a boundary CFT correlation 
function with several "probes" of the harmonic measure: 




-<\w'(z 1 )\ h K..\w'(z n+1 )\ h ^ 



n+l 



C = {WOy^Z^n+^ioo)) 



i=l 



This correlation function is found as 



4ai,n+ia: 



nl 1 4a'- a' 

\Zi-Zj\ * 1 



IK 




(60) 



When all Zi are at the same distance r from the origin it scales as 



(61) 
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As before, this correlation function is equal to the statistical average of a 
certain correlation function in the fluctuating domain, and we further apply 
the uniformizing map w(z) to transform this domain into the upper half 
plane H: 

i i 

Unlike eq. (|58p . the correlation function inside the average scales as: 
H(w( Zl )-^)f^ J] K^)-^)| 4q X- J] ( w ( Zl )-^-)) 2 < a '>, (63) 

i i<j i^j 



where a\ = = ao + y ag + K- We specifically separated the i = j terms 
since only they contribute to the short-distance behavior. All the other 
terms insure that the realizations of the curves in which any two points Z{ 
end up in the same sector are suppressed since the distances w{z{) — w(zj) 
are then small. Then the short-distance dependence of eq. (|62p is 

i 

Insofar as the scaling with r is concerned, we further approximate w(z{) — 
w(zi) ~ |zj| ~ r\w'(zi)\. This gives 

C ^E^JJ \ w >( Zi )fK+H V . (64) 
i 

We denote the exponents inside the average 

hi = 2h'i + 2a'*. (65) 

These are the weights of operators whose conformal charges = ao + 
\J a\ + hi are just 2a\. 

Finally, comparing eqs. ([6T]) and (f64"|) . we get the result 

<\w\z x )\ hl . . . \w'(z n+ x)\ hn+1 y oc r ^ (n) (hi,...h n+ i)^ 

with the higher multifractal exponent 

n+l ra+1 

A (n) (/ii, ...h n+ i) = 2q 1;n+ i ^ a hi + 2 ^ a hi a hj , 

i=l i<j 
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which is the formula (|55D . 

If Zi's are kept independent the computation becomes more complicated. 
It can be done to the end in the case of a single curve. In this case the points 
can be taken on the boundary and C reduces to a four-point function which 
contains tpi^' 

C(x 1 ,x 2 ) ^^w'{x l pw'{x 2 py ~ (O ftl (a:i)C» ft2 (x 2 )^ lj2 (0)*(oo))H, (66) 

where x± and x 2 are two points on the real axis. Since ^1,2 is degenerate on 
level 2, C satisfies the equation (|28p. Same result can be obtained within 
the SLE approach. Up to a normalization the solution is [3] 

C ~ xt^xf {h2 \ Xl - x 2 ) A - A ^)- A ^)Ff^^V (67) 

\x± — X 2 J 

A = /loo - ^1,2 - hi - h 2 , 

where /loo is the weight of VP(oo) and F{z) satisfies a hypergeometric equa- 
tion: 



K 

T 



z(l - + 1 + -A(M- 2 + -A(/ il ) + -A(/ i2 ) z 



2 v ' V 2 v ' 2 



- (A(fci,/i 2 )- A)F(z) =0. 

One of the two possible solutions of this equation diverges at small z and the 
other one is regular. If x 2 — > 00 the dependence of the correlation function 
on x\ should be C(xi,oo) ~ (O^ipi^ioo)) ~ xf^ hl \ We therefore should 
pick the regular solution of the hypergeometric equation. Then for C to be 
symmetric with respect to (xi,hx) <-> (x 2 ,h 2 ), F{z) should be a constant. 
A constant is a solution of the hypergeometric equation only if the weight 
at infinity is fixed by setting 

A = A(hi,h 2 ). 



We then obtain the formula (15 TH . 

We note that the function C(x\,x 2 ) as defined in eq. (|56"j) is real and 
positive for any x\, x 2 . Therefore, the constant F must have such a phase 
as to make the expression ([6T|) positive. 



6.3.4 Calculation of bulk multifractal exponents 

Calculation of the bulk multifractal behavior is done in much the same way 
as on the boundary, so we go straight to the general case of higher bulk 
exponents. The method of SLE is not yet developed for this case. 
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Let the critical system occupy the whole complex plane and be condi- 
tioned to have n critical curves growing from a single point in which we 
place the origin z = 0. We will assume that z = is the only common point 
of these curves, since the local results around this point are unaffected by 
the curves' behavior at large distances. The curves then form a star. We 
define the conformal map w(z) of the slit domain formed by the star to the 
exterior of a unit circle (the second picture in Fig. [TT|) and normalize w{z) 
at infinity as before: w(oo) = oo, w'(oo) = p~ 1 , where p is the conformal 
radius of the star. 

Close to the origin the curves divide the plane into n sectors. We consider 
a quantity 

■<\w'(z 1 )\ h K..\w'(z n )\ h ^, 

where z^s are points close to the origin, no two of them lying in one sector. 
As before, this is not a serious constraint because if some two points Z{ and 
Zj happen to be in the same sector, w{zi) — w(zj) — > when the size of the 
star L — > oo. 

Since n curves starting from the origin in the bulk are produced by the 
operator ^0^/2(0), we consider the correlation function 

n 

Cbulk = < JJO ft /(^,^ ,n/2(0)*(cX))>, 

i=l 

where h\ are related to hi as in f)65[) : hi have twice larger holomorphic 
charges: 

2a'i = a h . = a + ^"o + h i- 
The correlation function is found to be 

Cbulk ~ J] N 4q °'" /2Q * II \ Zi - z i\ A<a ' 3 ~ r 4Q0 -/ 2 ^^ +4 ^ a X', 

i i<j 

assuming that all Z{ are of the order of r. Proceeding exactly as in the 
boundary case we use the two-step averagin to rewrite the correlation func- 
tion Cb u ik through the conformal map of the fluctuating star: 

Cb U i k =^ni w '^)i 2h Kn°M(^^)'^)) xi/ ( oo )>c\s t a^ • 

i i 
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The correlation function inside the average is found to be the same as (j63[) 
although w(z) is defined differently. Combining the results we obtain 

-<|w / (^i)| /l1 . . . \w'(z n )\ hn y ~ r A ^ hl '- hn \ 
with the higher bulk exponent 

n n 

A b l(^i, • • • K) = + y, , 

i=l i<j 

where 

A bulk( /l ) = 2Q 0,n/2«/i - \f?h = ( 2 «0,n/2 ~ Ol Q )a h - ^ 

is the scaling exponent of a single <\w' (z)\ h >- in the presence of n critical 
curves in the bulk. These are the results quoted in eqs. ([541 [56]) . 

7 Conclusions 

In this work we attempted to clarify the subtle relation between the Gaussian 
Bose field description of critical behavior of statistical models and stochas- 
tic geometry of critical curves occurring in these systems. We further stud- 
ied the relation between two complementary approaches: the algebraic ap- 
proach of conformal field theory and a more direct approach of stochastic 
(Schramm-) Loewner evolution. 

A natural extension of this work is a study of global properties of critical 
curves such as the distribution of the area and harmonic moments of a closed 
critical loop. 

The critical curves considered in this paper all appear at 2D critical 
points described by CFTs with central charges c ^ 1. There are other 
CFTs which in addition to conformal symmetry have other symmetries, e.g. 
current algebras. These theories allow central charges c > 1. 

We have recently extended the SLE approach to the Wess-Zumino model 
with the current algebra SU(2)fc [39]. In this case the SLE trace carries an 
additional spin degree of freedom. The analog of the Gaussian field takes 
values in the Lie group. The critical curves which appear in such a system 
can also be studied by the methods described here. 
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